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Abstract 

We consider a quantum graph consisting of a ring with Rashba 
hamiltonian and an arbitrary number of semi-infinite wires attached. 
We describe the scattering matrix for this system and investigate spin 
filtering for a three terminal device. 



> 

(N ; 1 Introduction 

O 

Much work has been done on spin related transport properties of nano- 
j — ' electronic devices resulting in interesting applications, for example the so 

called spin field effect transistor proposed by Datta and Das [5]. There has 
c| \ been particular interest in using the Rashba effect to manipulate the spin 

de gree of freedom in such systems [H HH H9j [201 E21 [23] . In this paper we 
model a simple system exhibiting the Rashba effect, viz. a ring with Rashba 

i * i 

. hamiltonian attached to an arbitrary number of 'free' wires, using so called 

solvable models [H [211 El 02] • This means that we approximate the sys- 
tem by a — one dimensional — graph on which is defined an appropriate self 
adjoint Schrodinger operator. The advantage of this approach is that, as 
the name suggests, it allows us to get explicit expressions for the scattering 
matrix, and hence for the transport properties of the system, in this case in 
terms of the Greens function of the ring and the boundary conditions at the 
vertices. 

Our particular interest in considering this model is to investgate the possi- 
bility of constructing a spin filter. Various approaches have been taken to 
filter spin: we mention [8] in which the authors construct a spin filter using a 
four terminal device with the Rashba effect as well as |24j where the authors 
achieve spin filtering using a two terminal device and a magnetic field. A 
third approach, discussed in [13} 114j. uses a three terminal device with the 
Rashba effect and to some extent was the motivation for this paper. 
It is known that a device with two terminals and time reversal symmetry 
cannot polarise spin currents [15] (the device in [M] does not have time 
reversal invariance due to the magnetic field). Nevertheless, Kiselev and 
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Kim [13} [Tl] show that a three terminal device with time reversal symmetry 
and a particular geometric symmetry can make an effective spin filter. We 
consider the same geometry as considered in [T3], viz. a ring with three 
wires and symmetry with respect to reflection across the line defined by 
the 'incoming' wire. Whereas Kiselev and Kim assume the Rashba effect is 
localised at the 'incoming' terminal in our model the Rashba hamiltonian is 
present uniformly on the whole ring. Kiselev and Kim use a sophisticated 
numerical model of the system to calculate transport properties while our 
model is of course solvable. 

We believe that the formalism of solvable models offers, in general, advan- 
tages over numerical studies in that it allows us to derive explicit expressions 
for scattering properties thereby identifying principal features of the system. 
Ideally, these may be used to help optimise the design (for instance for spin 
filtering). In particular, for the three terminal device described above we 
investigate how the polarisation is related to the resonant eigenvalues on 
the ring, the Rashba coefficient and the angle of attachment of the wires. 
We observe, as did Kiselev and Kim, that this system may be used as an 
efficient spin filter. 



2 Quantum graph with Rashba hamiltonian 

We consider a ring shaped quantum waveguide where the width of the waveg- 
uide and the incident electron energy is such that the ring may be considered 
one-dimensional. Furthermore, we assume that there is structural inversion 
asymmetry [26] so that a Rashba term appears in the hamiltonian on the 
ring. Normalising the radius to one it can be shown [W\ [2"3~] that the hamil- 
tonian has the form 

where 

Id a 

D ° = — 7 Zm + ^ ar i 
i dd 2 

Or = &x cos(#) + o y sin(0) , 

9 G [0,2-71") is the local coordinate on the ring; a x , a y , a z , id denote the 
Pauli spin matrices and the unit matrix respectively; and a describes the 
strength of the Rashba spin-orbit coupling. The solutions of the eigenequa- 
tion, Hof = k 2 f, are 

/±,o(0, k) = e - ia * 9/2 e - ia ^ 12 e ±^±(^) (i) 

where k± = ^ k 2 + f ± \J\ + ^f and tan(^) = a, <p G (-f,|). The 
eigenvalues on the ring 

A±,„ = n 2 - (i ± n) (Vl + a 2 - l) 

correspond to the zeroes of cos(k±7t). Each eigenvalue A± jn has multiplicity 
two, the corresponding eigenspace is spanned by {e ±m ®x v e =Fme x i } where 

/ cos(^/2) \ / \ 

X t V e w sm(v/2) J ' x i \ cos(y>/2) J ' 
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Since A +in = A-,-^ and A +i „ < A_ >n we assume n £ {1,2,...} for X +jU 
and n £ {0, 1, . . .} for A_ >n . Finally, we note that the twofold degeneracy 
of the eigenvalues drops to a fourfold degeneracy when y/l + a 2 — 1 = m € 
{0, 1, . . .}. In this case we see that A_ 5 „ = A +i „ +m . 

Mostly we will write eigenfunctions with both spin eigenstates together in a 
2x2 matrix in order to simplify notation. In particular the solutions /± ; o 
may be used to find the Greens function, ie. the continuous solution of 



H a(e, m k 2 ) 



G(8, V ;k 



k 2 G(e,r];k 2 ) 
id 

G*(r],9;k 2 ), k £~R\a (H a ) , 



which is in fact 



G(o, m k 2 ) 

e -ia z K + T) 



(2) 



COs(K + 7r) 



f-,o(0) 



,-ia y (f/2 



cos(k_7t)J 2i(n + + k_) 



e ~ia z 9/2 e -i<r y <p/2 

2i(n + + k_) 



ia z K+{d—r\—-K) e —ia z K-(9—ri—n) 



COs(K + 7r) 



COs(k_7t) 



-ia y (f/2 ia z r)/2 



0", 



Here we take 6 — rj £ [0, 27r). 

We assume that the ring is attached to n semi-infinite wires. On each wire 
we have a 'free' hamiltonian 



H jfj - D j fj ' 
with generalised eigenfunctions 



1J_ 

i dx-i 



f ±? \dkxj 

J±,j — e 

where j £ {1, . . . ,n} is the index for the wire and Xj is the coordinate on 
the respective wire. 

We write the hamiltonian on the whole system 

n 

H = H ®Y J H J 

3=1 

and consider this as an operator on the Hilbert space L2(T, C 2 ) = L2(T, C 2 )© 
Sj=i ^2(1^+, C 2 ) of spinor valued functions on the graph T consisting of the 
ring T with n wires M + attached. To define this as a self adjoint operator 
we need to correctly define the domain of H which is related to self adjoint 
boundary conditions arising from the vanishing of the boundary form 



(Hf,g)-(f,Hg) = i^2({D j f,g) + {f,D j g))\ x . =Q 

j'=i 

+ iit(( D °f>9) + (f,D 09 ))\l = = f- ■ (3) 
3=1 1 
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Generally these boundary conditions are parameterised by a unitary matrix, 
for details see [H [3TJ. Here (•,•) is the inner product on L2(T,C 2 ), (•,•) 
is the inner product on spinors and {#j}™ =1 are the points where the wires 
are attached to the ring. 

We always assume that each vertex has three incident edges, ie. no two 
wires are attached at the same point on the ring. For a given vertex we 
denote by {*0i}? = i the values of the eigenspinor on edge i in the limit as we 
approach the vertex and by {VOi=i * ne va l ues of the outward derivative on 
edge i in the limit as we approach the vertex. In this paper we assume the 
following boundary conditions at the vertices 

p- = V2 = > M + V4 + V4 = o, (4) 

motivated by the fact that they are closely related to the ansatz for the 
scattering matrix of the T-junction as described in the physics literature [4] 
(we describe this relationship in the first appendix). Here edge % = 1 is the 
semi- infinite wire while edges i = 2, 3 are on the ring. The coefficient (5 
describes the strength of the coupling between the wire and the ring. 
These boundary conditions (jl]) are assumed to hold with the same (3 for each 
component of the spinor, ie. the coupling is independent of spin. However, 
j3 may in general be different at each vertex or point where a wire is attached 
to the ring. It is clear that these boundary conditions are self-adjoint, ie. 
([3]) vanishes (see [T7] for a discussion of boundary conditions in the presence 
of magnetic terms). 



3 Scattered waves and the scattering matrix 

The scattered waves tpi are eigenfunctions on the quantum graph satisfying 
the boundary conditions at the vertices and having the following form on 
the wires: 

4>i = f+,i + f-,iSu © ^2f-,jSji. (5) 

We reiterate that tpi corresponds to two spinor valued waves, one with spin 
up and one with spin down incident waves. Similarly, the components of the 
scattering matrix Sji , i, j £ {1, . . . , n}, are 2x2 matrix valued. Due to the 
nature of the boundary conditions we can define the scattered wave on the 
ring as 

^ i = Y / G k A ki (6) 

k 

where G k = G(9,9k), {9i}™ =1 are the points where the wires are attached 
to the ring and A is a matrix of coefficients. 

The boundary conditions for the scattered waves can be neatly expressed 
using the Greens function. Defining 

Gjk = G (Oj, 9 k ; k 2 ) 

we see from ([5]) that on the wires 

i>i\j = $ji + Sji , ip'ilj = ik (5ji - Sji) 
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where -| ■ denotes evaluation at the point where the j-th wire is attached 
to the ring and 5u should be interpreted as matrix with matrix valued 
components. Likewise from ^ we have 

k 

on the ring where we have used the property of the derivative of the Greens 
function. It is then easy to see that the boundary conditions in @ can be 
written 

I + S = PGA, ik/3 (I — S) = —A 

respectively where (5 is a diagonal matrix containing the coupling strengths 
for the n vertices. Solving for the scattering matrix we get 

s = (ikpgp + 1) (ikpgp - ly 1 . (7) 

This result can be generalised, at least in principle, using the techniques 
described in [10\ [T6] to find the scattering matrix of a quantum graph con- 
sisting of n semi-infinite wires attached to a compact graph consisting of m 
rings with Rashba term connected by edges of finite length. 



4 Spin filtering using a three terminal Rashba ring 

Let us consider a a three terminal device with symmetry as illustrated in 
figure [H ie. the angle £ G (0, 7r) between the first and second and first and 
third wires the same. To be precise we also need that the coupling constants 




Figure 1: The three terminal Rashba ring. 

at vertices two and three are the same — in fact, for simplicity, we will set 
all of the coupling constants equal to one, 0i = l. 

As was shown in |14j a three terminal device with this symmetry can poten- 
tially act as a spin filter. Specifically, for unpolarised current entering the 
first wire the polarisation of flux measured on wires two, P2i, a , and three, 
P31 aj along the a-axis satisfies 

-^21,1 = —P3l,x j P21,y = P31,y j ^21,2 = — -^31,2 • 
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A proof of these statements, following [H] but in the context of quantum 
graphs, is given in appendix 2. 

Appendix 3 contain an outline of the derivation of expressions for the con- 
ductance T21 and the polarisation in the z-axis P21.Z for current going from 
wire one to two for the device in figure [TJ In figures [2rEl we plot the conduc- 
tance T21 (upper curve) and Pi\, z (lower curve) against the energy A = k 2 . 
In these plots we assume a = 0.8; there is no significant change in the form 
of the T21 and P21,* curves with respect to a apart from at the special val- 
ues discussed below. We assume that the angle of attachment of the wires, 
£ = pir/q, is an integer fraction of tt. The resonance eigenvalues of the ring 



are indicated by □ for A +jn and for A_ jn (the leftmost □ and are at A +i i 
and A_ 5 i respectively with A_ 5 o < off the left of the plot). 
It is clear from these plots that at some of the resonant energies, depending 
on the value of £, this device would function as an efficient spin filter. Our 
case £ = 7r/2, figure most closely resembles the device considered in [H] 
and we observe as did Kiselev and Kim that the polarisation is maximum 
at odd resonant eigenvalues and severely damped at even eigenvalues — an 
explanation in terms of interference effects is given in [J3]. However, for 
our device at £ = tt/2 the polarisation changes sign at the odd resonances 
indicating that this may not be an ideal parameter regime for spin filtering. 
Clearly defined and isolated peaks can be found for instance at £ = 27r/3 
and we believe would be a better parameter regime for spin filtering. 
We remark that for £ = pir/q, T21 and Pi\ )Z ar e periodic repeating after 2q 
resonances A± jn with zeroes in T21 and Pi\ iZ at ^±,kq- It seems plausible that 
this behaviour can again be explained by spin interference and topological 
phase effects [15] . 

Finally, we observe from the expression for the polarisation P21 zi equation 
(|13p . the following behaviour as a function of the Rashba coefficient a. As 
noted above, for yl + a 2 — 1 = m G {0, !,•••} the eigenvalues on the ring 
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Figure 2: T21 and Pn,z for a = 0.8 and £ = tt/2. 
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become fourfold degenerate. For these values of a, k + = k_ + m + 1 so 
that — see the expression for Q given after equation (fT3j) — Q = and the 
polarisation is identically zero. This is clearly a generalisation of the degen- 
erate behaviour at a = 0, a possible physical explanation involves noting 
that there is an integral number of effective flux quanta [23] through the 
ring for these values of a and showing that for integral quanta the spin 
interaction is cancelled. We rather consider the gauge transformation 



on the ring where V{6) = e -*M™+l)0/2 e «W2 e ia z 6/2_ Reie ^ gince ig 
not single valued, we take 9 £ [0, 2tt) and assume that the point of at- 
tachment 6i = 0. We also make a change of basis on the wires, formally 
V(9j)HjV*{6j). Generally this results in all spin interaction being concen- 
trated at the boundary conditions; however, for Vl + a 2 — 1 = m G {0, 1, . . .} 
the boundary conditions remain spin independent so that, in analogy with 
a = 0, we have vanishing of the polarisation. Precisely, at all vertices except 
#1 = the original boundary conditions @ continue to hold while at the 
origin we get the new boundary conditions 



Here edge i = 1 is the semi-infinite wire, edge i = 2 corresponds to 9 G 
(O,^)) edge i = 3 corresponds to 9 £ (9 n ,2ir) on the ring and we recall 
that the above derivatives are in the outward direction from the vertex. 
These boundary conditions show that, up to an energy shift on the ring, the 
scattering properties of the ring fall into two classes depending on whether 
m is even or odd indicating that the conductance is a periodic function of 
a (as observed in the case of a two terminal device in |20j). 

Acknowledgements 

The author has benefitted greatly from conversations with Prof B.Pavlov 
and Dr U. Ziilicke. 

Appendix: Scattering matrix for the T-junction 

There is a well established description of the scattering matrix for the T- 
junction in the physics literature (see 0], pg 173, and [23 In this ap- 
pendix we show how this ansatz is related to the solvable models approach 
of specifying boundary conditions at the vertex of the T-junction. 
We note that (JH) are of 'projection type' [9], ie. we can express these bound- 
ary conditions in the form 




/rVi = ^2 = (-in 



- 3 , M + i>2 + (-irv 3 = o. 



P^iP = , Pip' = 
where $ = (V>i, tfa, ^) T , $ = (ip' 1 ,ip 2 ,tp' 3 ) T and 

1 / (3 2 (5 




(8) 
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P = I — P are projections. Now we suppose that the T-junction is at 
the point where three semi-infinite wires meet, instead of at the point of 
connection of one wire and a ring. As above we construct scattered waves 
for this non compact system and derive the scattering matrix which turns 
out to be [S] 



This scattering matrix is the same as posited in [H 1221 [25] for the T-junction 
(a different ordering of the edges is used in the last two references) motivating 
our choice of ([!]) for the boundary conditions. 

We note that, at least for the chosen boundary conditions ([HE]), scattering 
at the idealised T-junction is independent of energy. It is not difficult to see 
that this is true for any projection type boundary condition and furthermore 
means that there are no discrete eigenvalues or resonances associated to 
the junction (this is an equivalence [9], there are no discrete eigenvalues 
or resonances associated to the junction iff the boundary conditions are of 
projection type). We also note that projection type boundary conditions, 
also referred to in the degree two case as Fiilop-Tsutsui point interactions, 
are important in the study of 'quantum chaotic' behavior in quantum graphs 
|12j due to their scale invariance. 

Given a two or three dimensional quantum network, the problem of deriving 
appropriate boundary conditions at the vertices of an approximating one 
dimensional quantum graph is an active area of research (see (2j [7J [TT] ) as 
is the 'inverse problem' of constructing sequences of graphs with regular 
potentials so that in the limit we observe a chosen boundary condition from 
the whole U{n) parameter space [6, 3j. 

Appendix: Symmetries of the scattering matrix 

Here we follow the argument of [141 [T5] , which describes the symmetries of 
the scattering matrix, using terms appropriate for quantum graphs. 
We think of T as a graph in the plane and suppose that 7 is a closed curve 
nowhere tangent to T. The wronskian is defined as 



where {x{\ = 7 n T and cx(xj) is the orientation of the ordered pair formed 
of the orientation of the local variable at Xi and the orientation of 7 at X{. 
The operator D is one of Dq or Dj depending on whether X{ is on the ring 
or the wires. 

We always assume that the wronskian acts on solutions of the eigenequation, 
Hf = k 2 f, Hg = k 2 g, from which it is easy to see that W 7 (/, g) is piecewise 
constant. Furthermore, we assume that 7 is large, in particular it encloses 
and has no intersections with the ring, in which case we drop the subscript 




p 2 - 2 2/3 2/3 
2/3 -/3 2 2 
2/3 2 -ft 




W, (f,g) = J2(-lT (xi) ((Df,g) + (f,Dg))\ 
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and write 

W(f,g)=J2(( D ifi9) + (f,D i g))\ Xi . 

i 

From the constancy of the wronskian we see that on the ring 

n , 

£(<A)/,0> + </,A>0>)C'-=O 

3=1 3 

while on the wires 

«A/,<7> + (f,D i9 ))\ Xi = ((Dif,g) + (f,D i9 ))\ Xi=0 

so that W (/, g) is actually equal to the boundary form ([3|). In particular, if 
/ and g are eigensolutions satisfying the boundary conditions, ie. such that 
the boundary form ([3]) vanishes, then 

W(f,g) = (Hf,g)-(f,Hg) = 0. 

(In fact since our boundary conditions are 'local' we will have W 7 (/, g) = 
for any 7, but we do not need this.) 

The wronskian allows us to identify symmetries of the scattering matrix. 
Consider the wronskian of two scattered waves: 

= W{M j ) = Y,{{D^ j ) + {^ i ,D^ j ))\ Xk=Q 

k 

= J2 -k ((5 lk - S* k ) {8 kj + S kj ) + (5 lk + S* k ) (5 kj - S kj )) , 

k 

we get immediately 

S*S = I. 

We note that in this case, since ipi is matrix valued, the wronskian W (ipi, tpj) 
is properly thought of as a 2 x 2 matrix. 

Further symmetries of the scattering matrix may be found from operators 
commuting with the hamiltonian. Here we are mainly interested in the case 
where the graph, and boundary conditions, are invariant with respect to a 
reflection in one of the coordinate axes R : y <-> —y. It is then clear that the 
hamiltonian will commute with 1Z = a y R and we have the vanishing of the 
wronskian 

= W (K^,^) = {(DK^j) + (Ripi, D^))\ Xk=0 

k 

= {{5 ik a y - R{Si k )o y ) {5 kj + S k j) - (S ik a y + R(S* k )a y ) (5 kj - S kj )) 

k 

or 

o y R(S*)o y S = I S = a y R(S)a y . (9) 
Here a y is block diagonal with a y on the diagonal. 

Before we apply this we need to define some notation. It is convenient for 
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us to decompose the components of the scattering matrix in terms of spin 
matrices 

Sij — ^ ^ ] &aSij,a ■ 

a 

In terms of this decomposition we can express the conductance Tij and the 
polarisation in the cy-axis Pij^a for waves going from wire j to wire i as 

Tij — 2 Msjjji) + I ,x I \sij,y\ \sij,z\ ^ (10) 
Pij,a — 4^y (Sij^Sij^a -f- Sjj' jQ _iSjj iQ ,_|_i) . (1 1) 

Now we consider the three terminal device illustrated in figure Q] which is 
clearly invariant with respect to R. Then ([9]) implies 



where R(sij) = Si'ji. In particular, the polarisation satisfies 

Pil,x = 49 (S21,1S21, X + S21, z S 2 l,y) = ~Pzi,x 

P2\,y = 49 (s21AS 2 l.y + S2l, x S2l, z ) = P^l,y 
P2l,z = 49 (S21,IS 2 1, Z + S21,yS 2 l,x) = "^31,2 • 

Finally we note that for the system under consideration the time reversal 
operator K, = a y K, where K is complex conjugation, along with C = a z L, 
where L reverses the sign of a or equivalently <p, both commute with the 
hamiltonian. Proceeding as above these can be used to find yet more sym- 
metries of the scattering matrix (see |15| for a further discussion). 



Appendix: Derivation of polarisation and conduc- 
tance 

Using equations ([21 [7]) we see that we can express the scattering matrix for 
the device illustrated in figure [T] as 

S - U* I :, b :-> II: /> :> ! U 



where 



b 


Z\ 


Zl ]( b 


Zl 


Zl 


Z\ 


b 


Z2 \ \ z l 


b 


Z2 


Zl 


Z2 


b J V Zt 


Z~2 


b 



JJ = e i<T yfl 2 



Zl = JK 




COs(K_|_7r) COs(k_7t) J 

b = zq + if3~ 2 = k (tan(K + 7r) + tan(K_7r)) + if3~ 2 
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k = —k/2(n+ + K-), j = io z and we have assumed that all of the coupling 
constants are equal, = (3. 

From the form of the scattering matrix we see that 

S 21 = e -j£/2 e -i<W2 ( Sl + i(JzSz ) e iay<p/2 

where, using Maple, 

si + js z = h —^- ( z i b ~ z\z 2 ) (12) 
and, due to the form of Zi, the determinant 

D = b 3 — bz 2 z 2 - 2bz\z\ + z\z 2 + z\z 2 
is a complex scalar. 

It is easy to show that the terms due to U in the expression for S21 make no 
contribution to the conductance, equation (flQ|) . 



T 21 =2[\ Sl \ 2 + \s z \ 2 

and introduce a multiplicative factor into the polarisation in the z-axis, 
equation (fTT|) , 

P 21>z = 2icos(tp) (sis z - sis 2 ) . 

Using (|12p we write Si and s z in terms of z\ and b (here again the form of z\ 
is important) which gives us 

T21 = j^m 2 + N 2 ) N 2 - \ (b + 5) + *?*.)) 

and 

8cos(93) . ,_ 2 2 _ , 
'* = \D\ 2 3 y lZ2 ~ 1 ^ ' 

We again use Maple to get explicit expressions and simplify — here it is im- 
portant to cancel common factors cos _2 (k+7t) cos~ 2 (k_7t) which appear in 
the numerator and denominator to aid simplification and avoid numerical 
instability. At this step we also put (3 = 1. The expression for the conduc- 
tance and polarisation are then 

T 21 (k,t,a) = x2 + y2 , P 2 i rZ {k,t,a) = x2 + y2 , (13) 



where 



X = - (4k 3 + 3k) sin(«; + + k_) 7r + 4K 3 sin(K + + k_) (2^ - 7r) 

— 8k 3 sin (k + + k_) (£ — 7r) 

Y = — (6K 2 + |) COS (k + + K-) 7T — \ COS (k_|_ — K_) 7T 

+ 2k 2 cos (k + + k_) (2£ — 7r) + 4k 2 COS («+ + K_) (£ — 7r) 
R = K 2 + 4k 4 + TjK 2 [cos (2k + 7t) + cos (2k__7t) — cos (k+ + K-) £ 

— cos (k + + K-) (£ — 27r) — cos (k + (£ — 2-k) + K_£) 
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— cos (k+C + K - (£ — 27r))] 

+ 2k 4 [cos (k + + «_)(£- 2vr) + cos («+ + k_) (3£ - 2vr)] 
+ 4k 4 [cos (k+ + K-) £ + COS (/C+ + K_) 2 (£ — 7r)] 
Q = k 3 [cos (2k_7t) — cos (2k + 7t) 

+ cos (k+2£ + K_2 (£ - 7r)) - cos (k+2 (C - 7r) + K_2£)] 
+ 2k 3 [cos (k+ (£ - 27T) + - cos (k + £, + K-(£- 27r))] . 
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Figure 3: T21 and P2i, z for a = 0.8 and £ = ir/3. 




Figure 4: T 21 and P 2 i,z for a = 0.8 and £ = 2vr/3. 
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Figure 5: T21 and P^i,z f° r a = 0.8 and ^ = 7r/4. 




Figure 6: T21 and P21.2 for a = 0.8 and £ = 3tt/4. 
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